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D ' Abstract. The theory of plasma physics offers a number of nontrivial examples of par- 

tial differential equations, which can be successfully treated with symmetry methods. We 
' propose three different examples which may illustrate the reciprocal advantage of this "in- 

teraction" between plasma physics and symmetry techniques. The examples include, in 
particular, the complete symmetry analysis of system of two PDE's, with the determination 
of some conditional and partial symmetries, the construction of group-invariant solutions, 
' and the symmetry classification of a nonlinear PDE. 
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i 1 Introduction 

The theory of plasma physics offers a number of highly nontrivial examples of partial differential 
equations, which can be successfully treated with symmetry methods. The three different exam- 
, pies we are going to consider may clearly emphasize the reciprocal advantage of this "interaction" 



between plasma physicists and symmetry theorists: for the latters, plasma theory offers the op- 
portunity of testing their methods in concrete problems; for the formers, symmetry techniques 
provide new solutions to their equations. 

Our first example actually deals with the Liouville equation. This equation admits well 
known properties and a rich literature is devoted to it; there are however various similar but 
not equivalent forms of this equation, as we shall briefly point out, which admit solutions with 
different properties, of course. After some comments on the general peculiarities of our equation 
and of its symmetry properties, we shall provide some new, physically relevant, solutions. 

The second example deals with a system of two coupled PDE's. This system appears to be 
very rich in various and nontrivial symmetry properties. We shall provide the algebra of its 
exact Lie point symmetries and examine some properties of this (infinite dimensional) algebra. 
The system will be also compared with the properties of similar equations coming from mag- 
netohydro dynamics. Next we shall study some invariant solutions of this system, and discuss 
their physical aspects. This system also admits some relevant "weak" (conditional and partial) 
symmetries, as we shall show, leading to other families of solutions. In addition, we shall dis- 
cuss some peculiarities of the problem, considering, in particular, the effect of "truncating" the 
system (putting some physical parameters equal to zero) on its symmetry properties. 

The final example is actually a "symmetry classification" problem. We will propose indeed 
a nonlinear PDE containing two arbitrary functions and examine how symmetry properties 
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depend on the choice of these functions, also with the help of the notion of the equivalence 
group. The conclusion is that, apart from very special cases, only dilation exact symmetries 
are admitted by some well defined classes of functions. In addition, one nontrivial conditional 
symmetry is admitted, leading to a particularly simple reduced ODE, which can be easily reduced 
to quadratures. 

In this paper we will be concerned only with Lie point symmetries, see |H EH EJ E3 EE , 
under usual and standard definitions and assumptions, as stated e.g. in [^J[7j. For a different 
approach to symmetries of differential equations, using differential form methods, see the recent 
review paper (S]. 

For all details concerning plasma physics, we refer to [§1E3 an d to the papers quoted therein. 



2 Plasma kinetic equilibria described by Liouville equation 

Under suitable approximations and in appropriate experimental conditions, the configuration 
of an isothermal plasma embedded in a stationary magnetic field is 2-dimensional (planar: say, 
in the plane x, y). It can be shown (see ^U] and ref. therein) that such a configuration is 
completely described by a function u = u(x,y), which is the component of the magnetic vector 
potential perpendicular to the x, y plane, and which satisfies the classical "elliptic" Liouville 
equation 

V 2 u + exp(2n) = 0. (1) 

Let us write, first of all, its most general 1-dimensional solution u = u\{x) which depends only 
on x: 

u\{x) = — In ^— cosh(crr + d)j , c, d = const, (2) 

with c = 1, d = 0, this is called in plasma theory the 1-dimensional Harris sheet pinch solution. 

As is well known (see e.g. |1 II 1121 1X3*] ). equation admits the infinite dimensional algebra 
of Lie point symmetries generated by the vector fields 

d d d 

X = £(x,y)— +r)(x,y)— + C(x,y)—, (3) 
ox ay ou 

where £(x,y), rj(x,y) are arbitrary harmonic conjugated functions, and C( x ,y) = —d^/dx. 

It can be remarked for the sake of completeness, that the "hyperbolic" Liouville equation, 
analogous to (0), 

d 2 u d 2 u /n . 

dx~ 2 ~ W + = ° 

admits similar but not identical symmetry properties: in particular, in the hyperbolic case the 
functions £, rj in © can be arbitrary functions (see ^JE] and Vol. 1, p. 204]). 
Putting, in the case of equation l(TJ). 

(j) = (j){z) = £ + irj with z = x + iy (4) 



the function 4>(z) is then an arbitrary holomorphic function of the complex variable z = x + iy, 
and the symmetry property stated above for equation (0) implies that, choosing such a func- 
tion 4>(z), one can construct continuous families of solutions starting from any given solution 
to Q. 
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On the other hand, it has been shown by Liouville |13l 1161 1171 118j that the most general 
solution to Q can be written in the form 

u(x,y) = In ( 2 fl, 2 ) i z = x-Hy, (5) 

Vi + bOOr/ 

where the "generating function" 7 = 7(2;) is another arbitrary holomorphic function, and 72(2) = 
dj/dz. For instance, the choice 7 = exp(z) produces the above mentioned 1-dimensional Harris 
solution. The slightly more general choice 7(2) = kex.p(kz) + k, with k = \/k 2 — 1, £ 1, 
> 1, gives the solution 

y) = — In ( cosh(A; x) + (n/k) cos(fc y)) 

which describes a magnetic field configuration with a chain of magnetic "islands" analogous to 
"cat's eyes" vortex chain in hydrodynamics. 

It can also be noted that, choosing 7(2) = z a (a > 0) in ([5]). one obtains radial solutions 
to (H 



u = — In 



2 2,2 

r = x + y 



which present a logarithmic singularity at the origin unless a = 1. More precisely, it can be 
shown the following result |18| I19j . 

The most general solution to which is globally defined and satisfies the "normalization 
condition" 

exp(2u) dxdy < 00 

is the solution generated by 7(2;) = c(z — Zq), c £ R, c / 0, i.e. a solution radially symmetric 
with respect to a point (xo,yo) £ R 2 , of the form (with k = 1/c) 



«(a;, y) = In 



2|fc| 



k 2 + {x- x ) 2 + (y - y ) 2 



which is known as the Bennet pinch solution. 

Putting 7(2) = exp (/3(z)) , one obtains from (JSJ) this other useful expression of the general 
solution to (J2) 



u(x, y) = — In 



cosh(Re/?(z)) 

Mz)\ 



This formula and the above equation Q are the two expressions commonly used in the appli- 
cations in plasma physics. It can be noted, incidentally, that if one considers the other equation 

V 2 u-exp(2n) = (6) 

(or V 2 it + exp(— 2u) = 0) slightly different from (^Q), several different representations are known 
of its general solutions [111 I13j , but not all admit an analogous expression for our equation Q . 
For instance, equation (JHJ admits the two different solutions depending only on x 

u^\x) = — In I sinhx| and u^ 2 \x) = — In | sin a; | 

whereas the solution given in © is the most general solution to Q depending on x. One can 
also show that the elliptic Liouville equation in the form © does not possess solutions valid in 
the entire plane, while equation does possess such solutions. 
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For a full discussion about this (not marginal) point and a comparison between the possible 
expressions representing the various solutions of these Liouville equations, see [20]. It can be 
also noted that equation @ has the same symmetry properties as equation lfT|). 

Another peculiarity shared by all these types of Liouville equations is the possibility of stu- 
dying and finding their solutions by means of Backlund transformations. Also in this approach, 
elliptic and hyperbolic equations actually behave a little bit differently: see e.g. ^HlESESl- We 
do not insist on this topic, which goes beyond the scope of the present paper. 

Coming back to our equation let us consider for instance the Lie symmetry determined 

by 4>{ z ) = z<2 i i- e - 

X=(x 2 - y 2 ) — + 2xy— - 2x—. 

ox ay ou 

Starting from the Harris solution, we then obtain the following family of solutions, describing 
a continuously deformed family of curved plasma sheet configurations 



u(x, y) = — In 



[l - 2Xx + A V) cosh — 



Ar 2 



2Xx + A 2 r 2 

where r 2 = x 2 + y 2 , A G M is the group parameter. The generating function 7 is now 

7 = 70, A) = exp 



1 - Xz 



in this way, the essential singularity of the generating function of Harris solution is shifted from 
infinity to x = 1/A. 

Even more in general, one can perform an arbitrary holomorphic transformation 

z — > z = ip(z) 

and consider then the new generating function 7(2) = 7(^(2:)); thanks to (J5J), one obtains in 
this way other solutions to Actually, it can be shown that any solution to Q can be 
obtained in this way; this means in particular that any solution can be transformed (locally) into 
the 1-dimensional solution (j2j by means of a suitable holomorphic transformation z — > ip(z). 

The close relationship between the generating function 7(2) of any solution of Q arid sym- 
metry properties, can be also emphasized by the following result ^ 



Proposition 1. Let u = u(x,y) be any solution to 7(21) its generating function and Xq 
the Lie vector field which leaves invariant this solution: then 

j z (z)(f>o(z) = *7(*0, ( 7 ) 
where (j)o( z ) is the holomorphic function which determines Xq. 

For instance, for the above solution (J2J), one has 7 = exp(cz + c') (with c, d G M), and then 
from 0o = i/c, i.e. Xq = d/dy, as expected. 

To conclude this section, let us point out, among many other solutions, the three following 
examples, for their special interest to plasma physics (see [TU]^. 

Choosing as generating function 

7(2?) = (l+pz) 1 ^, p>0 

one obtains a one-parameter family of magnetic configurations with circular field lines, which 
includes the Bennet pinch for p = 1, ring-like structures for < p < 1, and which gives the 
Harris solution at the limit p — > 0. Instead, the solution generated by 



j(z) = pEri(z) 
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produces a sort of bar-like structure in the plasma current distribution. Finally, starting from 



z+p 

z) = exp z 

z — p 



one gets a solution which can be used to model some important features of the Earth magnetotail. 



3 A system of PDE's for magnetized plasmas 

The case we are now considering, which is the most elaborate and richest in different symmetry 
properties, is a system of two PDE's for the two functions u = u(x,y,t), v = v(x,y,t). 
The system is the following 

d 

— (u - V 2 u) - {u- V 2 u,v} + {u, V 2 v} = 0, 
J^(VV) - {u,V 2 u} + {v,V 2 v} =0, 

where the {•, •} is defined by 

{/ g } = ^L^l _ ^9_^]_ 

dx dy dx dy ' 

or, written in a more "symmetric" and compact form, 
d 

— (u - V 2 u ± V 2 v) + { v ± u, u - V 2 u ± V 2 v } = 0. (8) 
at 

In this system of PDE's the functions u(x, y, t), v(x, y, t) describe respectively the time dependent 
magnetic and electric potentials of a 2-dimensional planar configuration of a fluid plasma, 
embedded in a strong magnetic field orthogonal to the x, y plane and with a shear magnetic 
field in this plane. 

3.1 The Lie algebra of symmetries 

With the help of some appropriate computer package, e.g. |23U24lP5] . it is possible to show the 
following result: 

Proposition 2. The system © admits the Lie algebra of symmetries generated by the following 
operators: 

d d d d 

Xi = -ST, X 2 =y—-x—, X 3 = —, (9) 
at ox ay ou 

d d x 2 + y 2 d . , 

x ^- Uj T, +tx T v + ^ai (10) 

and by the family of operators, depending on three arbitrary smooth functions H(t), A(t), B(t): 
X H = H{t)^ (11) 



X {A , B) = A(t) A + B (t)|. + - y At )^-, (12) 



where A t = dA/dt, etc., which generate an infinite dimensional subalgebra (actually: an ideal). 
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The nonvanishing commutation rules between the above defined symmetry generators are 
the following: 

[X 1 ,X 4 ] = -X A , [X 1 ,X H ] = X Ht , [Xi,X(a jB )] = X(A t ,B t )i 

[X 2 ,X( AjB )] = X(_ BjA }, [X4, X^b)] = — -X"(_ tB) tA)i [X(A,B),X(c,D)] = X H' 

where H = AC t - BD t - CB t + DA t . 

The operators X\, X2, X% generate respectively time translation, space rotation and the 
translation u — > u + k, and together with the operator X4 form a subalgebra of the algebra of the 
symmetries; the operators X^a^b) include in particular space translations d/dx and d/dy. The 
operator X B generates the translation v — ► v + H(t); this operator and X3 (which is the center 
of the algebra) do not change either the magnetic field or the plasma velocity. Some properties 
of the two remaining operators X4 and X^ B ) wm be considered in the next subsection. 

3.2 Symmetry properties and group-invariant solutions 

Before considering group- invariant solutions of system (jHJ) , let us look for an optimal system (see 
e.g. j2J) of 1-dimensional subalgebras of the algebra of symmetries listed in Proposition 2. It is 
easily seen that in the adjoint representation the action of the operators Xjj and m on all 
the vector fields X a (where X a denotes any of the operators defined above) has the effect of 
transforming the functions H(t), A(t), B(t) and does not touch the other operators; the same 
happens for the action of any X a on X B and X^ B ) ■ For instance, one has: 

Ad(e X p(AX (A ,B)))^ = ^ + AX (Ati B () - ^{AA tt - BB tt )^-, 

Ad( exp(AX (j4i B)))X (CiD) = X (A)B) - X(AC t - BD t - CB t + DA t )^-. 

In view of the physical applications, it is more convenient to consider the functions H, A, B 
as arbitrary functions, and therefore we do not include in the classification the operators X B , 
X(a,b)- Restricting to the operators X\, X2 and X4 (notice that X3 acts trivially on all the X a ), 
an optimal system is 

d d d d 

Xi = — , X 2 =y- x—, X a = a— + X 4: , 

at ox oy at 

where a is a real parameter. 

The equations for time-invariant or rotationally-invariant solutions can be easily deduced; on 
the other hand, they are actually very special solutions, of no general interest. The same is true 
for solutions invariant under X B (he. solutions with v = 0). We then consider in some detail 
the symmetries X4 and X( A ,b) ■ It i s convenient to examine first the symmetry Xrj^m . 

1) The symmetry X( A ,b) expresses the property that, if u(x, y, t), v(x, y, t) is a solution of (JSJ, 
then also 1 

u(x, y, t) := u(x - A(t),y - B(t),t), 

v(x, y, t) := xB t - yA t - -(AB t - A t B) + v(x- A(t), y - B(t),t) 

solve our system (jHJ), for any A(t), B(t). It must be noticed that this corresponds to a change 
of spatial coordinates into a moving frame which produces in turns the additional term xBt — 



x The parameter A which should be introduced to parametrize this family of solutions can be clearly absorbed 
in the (arbitrary) functions A and B. 
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yAt — (1/2) (ABi — AtB) in the component v. Since v is proportional to the electric potential, 
we can interpret this symmetry as expressing the fact that a time dependent, spatially uniform, 
electric field imposed on the system induces a uniform time-dependent electric drift. 

We want now to look for the solutions to (JSJ) which are invariant under Xi^m . Let us recall 
that if a system of PDE's for the q functions u a = u a (x) of the p variables Xi admits a symmetry 
X = iid/dxi + (ad/dua, i = 1, . • • , p, a = 1, . . . , q, then the X-invariant solutions must satisfy 
the q equations 

c„-e,^ = o. (i 3) 

These equations in our case take the form 

Au x + Bu y = 0, Av x + Bv y = xB t - yA t . (14) 

The first equation can be interpreted as the requirement that the displacement produced by the 
electric drift is parallel to the field lines of the shear magnetic field, while the second expresses 
the requirement that the electric potential of the displaced plasma element remains constant. 

Instead of looking directly for the invariant solutions to (JEJ starting from (|14j) , it can be more 
interesting to give preliminarily the expression of equations (jSJ) when written in the "canonical 
coordinates" (or "symmetry-adapted coordinates" ) [2J [7] . In our case these coordinates can be 
chosen as follows 

A(t)x + B(t)y 

° = B(t)*-A(t)y, w = -^_W|, 

where s is (together with t) an invariant coordinate under the symmetry X^ t ^, and w is the 
coordinate "along the action" of the symmetry, i.e. X^ ^w = 1. Notice also that {f,g} x ,y = 
{/^In- 
putting, as indicated by (|T4"|) . 

u = U(s,w,t), v = Q(x,y,t) + V(s,w,t) (15) 

where Q is a quadrupolar term given by 

Q = A 2 + B 2 (( A t B + AB t) i x2 ~ V 2 ) ~ 2xy(AA t - BB t j) , 
we obtain from (jHJ) the following system 

^-(U -V 2 U ±V 2 V) + \V ±U,U - V 2 U ±V 2 V\ 
ot 

+ TTo^r^^-iU ~V 2 U ±V 2 V) =0, (16) 
(A 2 + B 2 ) 2 dw y ' 

where V 2 = (A 2 + B 2 )d 2 /ds 2 + (A 2 + B 2 )~ l d 2 jdw 2 . Here the different role of the two variables s 
and w is evident: in particular - as expected - there is no explicit dependence on w, whereas 
the last term, proportional to s, can be interpreted as a Coriolis-type contribution. 
Now, if one looks for invariant solutions, i.e. for solutions of the form 

u = U (s,t), v = Q(x,y,t) + V (s,t) 

the above equations (|16|) become 



1(U - (A 2 + B 2 )U , SS ) =0, | ((A + B 2 )V 0ySS ) = 0. 
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These clearly imply 

Uq - {A 2 + B 2 )U 0>SS = F(s), {A 2 + B 2 )V 0>SS = G(s), 

where F(s), G(s) are arbitrary functions. Notice that d 2 /dx 2 + d 2 /dy 2 = (A 2 + B 2 )d 2 /ds 2 
and that these equations are actually ODE's, indeed the variable t here appears merely as 
a parameter. 

Special physically relevant solutions of this reduced system (and of the system (JHJ), of course) 
can be easily obtained starting from particular choices of the arbitrary functions A(t), B(t), 
F(s), G(s) introduced above (see [§] for some explicit examples). 

2) We now consider the symmetry X 4 of our system (JHJ): it implies that if u(x,y,t), v(x,y,t) 
is a solution, then also 



u{x,y,t) := u(x cos(At) + ysin(Ai), — xsin(At) + ycos(\t),t), 
v{x,y,t) := v(xcos(Ai) + ysin(At), — xsin(Ai) + ycos(At),i) + A 



x 2 + y 2 



is a family of solutions to (jSJ) for any A € M. This represents a sort of rotated solutions with 
angular velocity A plus a radial term in the component v which gives the additional velocity 
field corresponding to the rotation. 

The solutions to (jSJ), which are invariant under symmetry must satisfy the invariance 
condition 

du du dv dv r 2 

^ dx X dy ' ^ dx % dy 2i' 

where r 2 = x 2 + y 2 . It is easy to find that they are of the form, with = arccos(x/r), 

2 

u = U^(r,t), v= r —e + V { ~°\r,t), (17) 

where U^°'{r,t), V^°'(r, t) satisfy the linear and uncoupled homogeneous PDE's 

r 2 U^. - 2rtU^ t + rU^ - 2tU^ - r 2 U^ + 2rtuj; 0) + 3U® = 0, 
2rtV}$ - r 2 V$} + 2tv£ ] - rV^ + = 0. 

A simple solution of the first equation is, e.g., = r 2 t, whereas the second one admits 

solutions of the form V^- ' = r a t h with b = (a 2 — 2a — 4) /2a, for any a. It can be noted in 
particular that the special form of the component v in (|17[). which has a fixed term containing 
a "cut" discontinuity, looks as a spiral and expresses the fact that this solution contains a 9- 
independent azimuthal electric field, vanishing for t — ► dboo. Solutions of this form are not new 
in plasma theory and in magnetohydro dynamics: see e.g. [^1 an d Vol. 1, p. 393]. See also 
Subsection 3.4 for other remarks and comments on various related physical aspects. 

It is remarkable that, in both cases of symmetries (|12[l and (|1U(I . the invariant solution u 
satisfies a linear equation (therefore linear superposition principle holds), and that the same is 
true for v, apart from a fixed additional term (for some comment on this point we refer to [HI ) - 

3) Finally, one can also look for solutions invariant under symmetries of the form, e.g., 

d d 

M + X H , m +X {AtB) (18) 

(or similar combinations: see our remark on the optimal system of symmetries at the beginning 
of this section). Clearly, the most convenient choice of symmetry to be considered, as well as 
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the choice of the arbitrary functions H, A, B involved, can be suggested by the specific case to 
be examined or by the experimental configuration. Just to give some examples, let us consider 
the two symmetries (|18|> . 

For the first symmetry, it can be easily seen that the invariant solutions have the form 

u = U{x,y), v = T(t) + V(x,y), 

where Ti = H (t) and U, V satisfy the system 

{U - V 2 U,V} = {17, V 2 V}, {U,V 2 U} = {V,V 2 V}. 

Considering instead the symmetry d/dt + Xug), one finds that the invariant solutions have 
the form 

u = U(x-a(t),y-fi(t)), 

v = ( x -a(t))B- (y - fi(t))A + aB - fiA + V((x - a(t),y - fi(t)) , 
where at = A, fit = B, and U, V must satisfy a system of nonlinear PDE's. 



3.3 "Weak" symmetries 

The above system (JSJ admits also several interesting "non-exact" (here generically called "weak") 
symmetries; we are going to consider an example of conditional symmetry (see [271 1281 1291 13UI 
I3T1 13*2] ) , and one of partial symmetry (according to the definition given in (SHj , see also [31] and 
below). 

As well known, any vector field X = ^d/dxi + Cad/du a is a conditional symmetry for 
a differential equation (or a system thereof) A = if the system A = enlarged with the 
invariance condition (|13|) admits some solution. If this is the case, a conditional symmetry 
allows the reduction of the initial equation into a reduced form, and in this way one can obtain 
other invariant solutions: see [""71 [2EJ \M ES3 ED E2] • See also 0H1 E3 EH 1^1 SOI SIl 1121 S3] 
for careful discussions about various related problems and reduction procedures. 

An example of conditional symmetry leading to interesting solutions to (JSJ) is the following 
(more precisely, it is a "contact conditional symmetry"): 



d d 
ou ov 



From this symmetry, one finds e.g. solutions to (JSJ of this form 

u = sin [k(x - T(t))] , v = exp(±ny) - T t (t)y, 

where fc,K£K with k 2 — k 2 = 1 and T(t) is arbitrary, or similar solutions where the functions 
sin and exp are exchanged. Solutions of this form are interesting due to the presence of terms 
depending on x — T(t) describing generalized wave propagation. 

It is known that there are some delicate points related to the definition of conditional symmet- 
ries: see [SSJ EH] EZJ EH EM EI3 El H2] ; this is actually related to the introduction of a subtler 
classification of the notion of conditional symmetry |441 145] . We do not deal here with this 
problem, and we prefer to consider a particularly interesting example of partial symmetry for 
our system 0. 

While exact symmetries transform any solution of the given problem into another solution, 
conditional symmetries do not transform - in general - solutions into solutions, as well known. 
Partial symmetries play in a sense an intermediate role, indeed they transform solutions belong- 
ing to a proper subset of solutions into solutions in the same subset. This subset is defined 
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in the following way. Let us assume that some vector field X is not an exact symmetry for 
a differential equation (or a system) A = 0: this means that, denoting by X* the appropriate 
prolongation of X, one has X*(A)|a=o 0- Let us then introduce the condition 

A^ := X*(A) = (19) 

as a new equation, and consider the enlarged system 

A = A (1) = 0. (20) 

It is clear that if X is an exact symmetry of this enlarged system, then the subset S^> of the 
simultaneous solutions of this system (if not empty, of course) is a "symmetric set of solutions" 
to A = 0, i.e. a proper subset of solutions which have the property of being mapped the one 
into another by the vector field X. It is also clear that this property is not shared by the other 
solutions to A = not belonging to the subset SW. In principle, this procedure can be iterated, 
see |331 134j , but no iteration is necessary in the example we are considering. 
The vector field we want to deal with is 

X - u^. (21) 

First of all, notice that this is trivially a conditional symmetry for (JHJ, indeed the invariance 
condition simply amounts in this case to look for the special solutions to (jHJ) with u = 0, i.e. to the 
hydrodynamic limit. More interestingly, we now show that the vector field 1)21 Jl is a nontrivial 
partial symmetry: notice that considering this vector field as a partial symmetry corresponds to 
looking for solutions to (jHJ with the property that the component u admits a scaling property, 
i.e. for solutions u, v such that also Au, v solve (jHj) for A £ R. Applying the prolongation X* to 
the system (JHJ, and combining the resulting equation A^ = X*(A) with Q, one gets the new 
condition 

{u,V 2 u} = 0, (22) 

which is then the condition characterizing the subset of solutions with the above specified 
property. It is easy to verify that the system of the three equations (jHJ and (|22j) is symmetric 
under (|21j) . showing that ()21|) is indeed a partial symmetry for (jHJ. Thanks to (J22J, we deduce 
this equation for v 

l(V 2 v)+{v,V 2 v} = 

which does not contain u, and is equivalent to the two-dimensional Euler equation for an in- 
compressible fluid. On the other hand, the equation for u is linear in u: obviously, this agrees 
with the presence of the (partial) symmetry given by (|21j) . 

Before considering examples of particular simple solutions, let us come back to an important 
physical feature of our initial system of differential equations (jHJ. This system indeed would 
actually contain some physical parameters that we have normalized to the unity up to now. 
There are physical situations, however, where one of these parameters turns out to be negligible 
and can be put equal to zero with a good approximation [§]: this coefficient multiplies the 
term {u, V 2 u} in (JHJ. For this reason, it can be interesting to repeat calculations determining 
symmetry properties and solutions of the approximate system. The quite surprising result is 
summarized in the following 



Proposition 3. The truncated system 

— (u- V 2 u) -{u- V 2 u,v} =0, — (V 2 v) - {u,V 2 u} + {v,V 2 v} = (23) 
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admits precisely the same exact symmetries as the original one (JSJ). The same is also true if 
one or both of the other equations 

{u,V 2 v} = 0, {u,V 2 u} = (24) 

are appended to the above system ( therefore, even if the partial symmetry ()21[) is taken in con- 
sideration also within this approximation). 

Now observing that condition (|22j) implies V 2 u = A(u, t) for some smooth function A, and 
similarly {u, V 2 v} = implies V 2 t> = B(u,t), the system (|23|) - (j24|) becomes 

(1 - A u ) (u t - {u, v}) = A t , V 2 u = A{u, t), 
-B u (u t -{u,v}) = B t , V 2 v = B(u,t). 

Therefore, assuming e.g. A u = 1 forces At = 0; instead if A u ^ 1, B u ^ 0, but At = Bt = 0, the 
above system takes the very simple form 

u t = {u,v}, V 2 u = A(u), V 2 v = B(u). (25) 

Elementary solutions of this system (which are also solutions to (JSJ), of course) can be imme- 
diately found, for instance, just to give some simple examples, 

u = ci sin(k(x - t)) + c 2 sm(k(y - t)) + c 3 ), v = x - y, 

where Cj, k are arbitrary constants, or also u = 2t — 6, v = x 2 + y 2 , with 6 = arctan(?//x), and 
u = U(y — t), v = x, where U is an arbitrary regular function. 

3.4 General remarks 

As a general remark about the problem examined in this section let us point out that the set 
of equations © are of the interest for the study of the nonlinear dynamics of fluid plasma 
configurations and in particular for the interaction of magnetized plasma vortices and for the 
development of collisionless magnetic field line reconnection. In particular, the explicit solutions 
found in the present article using Lie point symmetries have a direct physical interpretation; for 
instance, the solutions l|15 |l — (|16 |) have a role in the recent investigation of nonlinear evolution and 
saturation of magnetic reconnection instabilities |46|l47j as plasma configurations forced from the 
boundaries. We note that equations (JHJ) obey boundary conditions that are determined by values 
of the electric and magnetic fields at the boundaries of the domain where the plasma is enclosed. 
These boundary conditions may include plasma fluxes through the boundaries. Referring for 
example to the solutions ()15|) - (|16|) . we see that they correspond to a plasma velocity pattern 
controlled by the boundary conditions that determine the arbitrary functions A(t) and B(t). 
The required boundary conditions correspond to imposing the electric field at these plates: 
the tangential component is proportional to Vv and determines the velocity pattern, while the 
normal component is proportional to du/dt and determines the magnetic terms. 

It is interesting to observe that the system (jHJ) can be seen as a generalization of the two- 
dimensional limit (i.e., the case of solutions independent of the coordinate z orthogonal to the 
plane (x,y)), of the so-called Reduced Magnetohydrodynamic equations (Kadomtsev-Pogutse 
equations) [IE], see also The symmetry properties of these equations, that are frequently 
used to describe plasmas embedded in a very strong magnetic field, have been studied in 
see also Vol. 2, p. 390]. It can be remarked that a common feature of the vector fields 
describing the symmetry properties of all these equations is the presence in their expression 
of one or more arbitrary functions (see also |49| ^ constant along magnetic (and velocity) field 
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lines. The symmetry properties of the three dimensional visco-resistive magnetohydrodynamic 
equations in the incompressible limit have been studied in [5U], see also (HI Vol. 1, p. 389]. 

Finally, the study of the partial symmetry (|21j) has been useful not only to obtain a reduction 
of the initial equations, leading also to new particular solutions, but also to point out a special 
property of the system (jHJ): it can be suitably truncated as in equation 1)23(1 (i.e. a physical 
parameter can be put equal to zero) without altering its general symmetry properties. 

4 Symmetry classification of an equation 
for axisymmetric plasma configurations 

Our third example deals with a completely different situation. The equation we are going 
to examine contains indeed two arbitrary functions F(u) and G(u) of the unknown variable 
u = u(x,y), and the goal is now to perform the symmetry classification of this equation, i.e. to 
find those F, G for which the equation admits nontrivial symmetries. In general, the symmetry 
properties of an equation may strongly depend on the choice of the arbitrary functions involved. 

Just to recall briefly a significant example, let us point out the case of the nonlinear Laplace 
equation V 2 u = F(u), with u = u(x,y) (or its hyperbolic counterpart, i.e. the nonlinear wave 
equation u xx — u yy = F(u)). Symmetry properties of equations of this form have been exensively 
studied; a complete symmetry classification of the nonlinear wave equations can be found in |14j . 
Summarizing the results, one has that these equations, excluding for a moment the particular 
cases F = u and F = 1, admit the infinite dimensional algebra of symmetries already mentioned 
(Section 2) if F(u) = exp(u) (the Liouville case), the nearly trivial symmetry 

X = (k-1) ( X ±+yjp\-2u£- 

\ ox Oy J ou 

if F(u) = u k , and only the obvious symmetries (translations and rotations of the variables x, y) 
otherwise. In the case F = u one finds the standard symmetries of the linear equations; finally, 
if F = 1, the infinite-dimensional algebra of symmetries given by ® is replaced by 

d d d 

x = v) dx~ + y) dy~ + ( cu + A ^ X: ^Thl' 

where c = const and A(x,y) must satisfy 
V 2 A - 2£ x + c = 0. 

Performing the symmetry classification of a given equation may be not an easy task: for some 
explicit examples and the presentation of the techniques used in each situation, see e.g. [^D E21 
E31 02J EH] and the papers quoted therein. 

The PDE we want to consider is 

Uxx + ^Ux + u yy = x 2p F(u) + G(u), (26) 

where a, p are constants. With a = — 1, p = 1 this equation is known in plasma physics as 
the Grad-Schliiter-Shafranov equation (see [HE]) and describes the magnetohydrodynamic force 
balance in a magnetically confined toroidal plasma. In this context, u is the so-called magnetic 
flux variable, x is a radial variable, then x > 0, while the two arbitrary flux functions F(u), 
G(u) are related to the plasma pressure and current density profiles. 

According to the standard definitions and procedure (see pQ), we can now look for the equi- 
valence group; preliminarily, we look for the kernel of the full groups of equation (|26|) . i.e. the 
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intersection of all groups admitted by Q26JI for any arbitrary choice of F and G. It turns out 
that this kernel is almost trivial: it contains indeed only the translation of the variable y. 

The complete symmetry analysis of the above equation ()26|) is a little bit tedious, but some 
simplification is provided for equations of this type, where the nonlinearity is present only in 
the arbitrary functions F, G. In this case indeed one can show that the coefficients £, 77, £ 
of the admitted symmetries X = £d/dx + r/d/dy + Qd/du are subjected to some precise linear 
constraints, which admit a nice geometrical interpretation j57| . 

In the case of our equation (|26[) . one obtains that, apart from the translation of the variable y, 
the equivalence group contains the following transformations: 

i) the translation of the dependent variable u, i.e. u — > u + k, 

ii) the scaling 

u — ► au, F — > aF, G — > aG, 
m) the scaling 

x^ra, y^ay, F^a" 2(p+1) F, G -»■ cT 2 G. 
The conclusion can be summarized as follows. 

Proposition 4. Apart from the particular cases where both F and G have the form c\ + c^u, to 
be examined later, the only choices for the functions F, G which lead to an equation admitting 
some symmetry (up to transition to equivalent functions via equivalence group) are the following: 

a) F(u) = u 1+ ( p+1 )^, G(u) = u l+1 / q 

for all q 7^ 0, where the admitted symmetry is 

d d n d 

X = x— + y- ^qu— 

ox oy au 

and 

b) F{u) = exp ((1 +p)u), G(u) = exp(u) 

with symmetry 

d d d 
X = x— + y— -2— . 

ox oy ou 

In addition, one has that 

c) if in equation (|2li|) p = — 1 and G(u) = 0, then the scaling of the variables x, y, generated 

by 

d d 

X = x 7T + y7T 
ox oy 

is a symmetry for all F(u) . 

For instance, with account of the equivalence group transformations, the most general choice 
of functions in class a), is 

a') F(u)= Cl (c + u) 1+{p+1)/q , G(u)=c 2 (c + u) 1+(1/q \ (27) 



where c, ci, C2 are constants, with symmetry 
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No new symmetry appears with particular choices of the parameters a and p, and of the 
functions F, G as well, apart from the following almost obvious exceptions. First of all, we 
exclude from our classification the case a = and p = (or a = and F(u) = 0) because in 
this case our equation coincides with the already mentioned and well known nonlinear Laplace 
equation. Other possibilities appear if both F and G have the form c\ + c 2 u. If F = G = u, 
the equation exhibits the standard symmetries of any linear equation. If F = 1, G = u, there is 
a new admitted symmetry which can be written in the following way: 

X = {u + ^(x,y))^-, (28) 

where, with S (u) denoting the left hand side of equation (|26j) . $ solves the equation 
£($) = $ - x 2p . 

If F = u, G = 1, the new symmetry has the same form as (|28|) . but <I> must solve 
£($) = x 2p § - 1. 

Thanks to the equivalence group of our equation, it remains only to consider the case F = G = 1. 
In this case, although of no great interest in the applications, the situation is much more com- 
plicated. Indeed, the most general vector field X belonging to the algebra of the admitted 
symmetries can be globally described in the following way: 

X = (cix + 2c 2 xy)-^ + (c + ciy - c 2 (x 2 - y 2 ))^ + ((c3 - c 2 ay)u + y)) J^, 

where the function fy(x,y) satisfies the equation 

+ c 3 - 2ci + x 2p (c 3 - 2 Cl (p + 1)) - c 2 y(4: + a) - c 2 x 2p y(4(p + 1) + a) = 0. 

There is, in addition, an interesting conditional symmetry for equation (|26|) in the case 
a = — p (which covers in particular the concrete case a = —p = — 1 of the Grad-Schluter- 
Shafranov equation), and G(u) = k 2 F(u) (k = const). The conditional symmetry is given by 

,9 „ d 
X = k— + x p — . 
ox ay 

This symmetry is particularly interesting and useful because, with introduction of the symmetry- 
invariant variable 

p+1 — ky if p — 1, or s = In x — ky if p = — 1 



p + 1 

the PDE is transformed into the very simple reduced ODE 

u ss = F(u) 

which can be easily integrated. 

The class (a) (or (a')) symmetry in Proposition 4 corresponds, for a = — p = —1, i.e. for 
the Grad-Schliiter-Shafranov equation, to spatial dilations where u scales as x~ 2q , y~ 2q . An 
elementary solution obeying such a symmetry is obtained by taking, e.g., q = —2 and by 
considering "cylindrical" solutions with du/dy = 0. Choosing F(u) = c\ and G(u) = (8— c^u 1 / 2 , 
according to equations (|27|). we obtain from ()26|) u = x 4 . By recalling the physical meaning of 
the functions F(u) = c\ and G(u) = (8 — ciju 1 / 2 we find that this solution describes a plasma 
configuration in the (dilation dependent) domain < x < xq with a pressure p of the form p = 
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(ci/47t)(xq — x 4 ) and a total current I along the cylinder of the form I 2 = [4(8 — ci)/3](x[j — x 6 ), 
where [4(8 — ci)/3] 1 ' 2 Xq is the current flowing along a conductor at x = 0. The boundary 
conditions have been chosen such that the pressure and the total current vanish at x = xq. 

Finally we remark that the Grad-Schluter-Shafranov equation is a special example of a family 
of equations that are used in plasma physics and in space physics and astrophysics [HH1 in 
order to describe time independent configurations of a plasma embedded in a magnetic field in 
the presence of plasma fluid motions (not included in the Grad-Schluter-Shafranov equation). 
Similarly to the Grad-Schliiter-Shafranov equation, these equations depend on a set of arbitrary 
functions, but can change from elliptic to hyperbolic type depending on the (local) value of the 
plasma velocity. 

5 Conclusion 

In conclusion, we believe that the three examples considered in this paper may give a detailed 
idea of the opportunities offered by symmetry methods to the study of differential equations 
arising from modern physical theories. Let us emphasize in particular the completely different 
peculiarities of our examples, all suggested by plasma physics, and correspondingly the different 
role played by symmetry techniques in their application to each situation. For other different 
applications to plasma physics of similar methods and procedures, based on general symmetry 
properties, we refer also, e.g., to the papers |59 j 160 } l6T| 162]. 
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